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1 Introduction 



In this paper, we compare recent classical and quantal derivations of induced gluon ra- 
diation for applications to ultra-relativistic nuclear collisions. The classical distribution, 
based on the McLerran-Venugopalan model [|I|, was recently computed to order in 
Ref. 0. The soft gluon bremsstrahlung distribution was computed via pQCD by Bertsch 
and Gunion in Ref.|^ within the Low-Nussinov approximation. Another quantal distri- 
bution, based on the Gribov-Levin-Ryzkin ladder approximation[^, was recently applied 
by Eskola et aip|. Finally, there has been considerable recent effort to compute mod- 
erate j9x (mini-jet) distributions based on the conventional collinear factorized pQCD 
approach [g % %. 

Interest in the moderate gluon distributions arises in connection with estimates 
of the initial conditions and early evolution of the quark-gluon plasma formed in ultra- 
relativistic nuclear collisions at RHIC (^s = 200 AGeV) and LHC (^s = 6500 AGeV) 
energies. Until recently, the main source of mid-rapidity gluons was assumed to be copious 
mini-jet production as predicted via the conventional pQCD gg gg processes!^, |^, ^. 
However, in Refs.[||, 0] it was suggested that another important source of mid-rapidity 
gluons could be the classical Yang-Mills bremsstrahlung associated with the passage of two 
heavy nuclei through each other. In the conventional approach, beam jet bremsstrahlung 
is assumed to influence only the non-perturbative low transverse momentum beam jet 
regions. Beam jets are then typically modeled by pair production in Lund or Dual Parton 
Model strings. See for example Refs. 0, |^ and references therein. 

The novel suggestion in was that for sufficiently large A nuclei and high energy, 
the initial nuclear parton density per unit area could become so high that the intrinsic 
transverse momentum of the partons, ^/x A^^^Aqcd, could extend into the mini-jet 
perturbative regime, /c^ ~ 2 — 4 GeV. It was suggested that beam jet bremsstrahlung 
could even dominate that few GeV transverse momentum region because it is formally of 
lower order in as than mini-jet production. Such a new source of moderate p± partons 
would then significantly modify the early r ~ evolution and hence possibly modify 

many of the proposed signatures of the quark-gluon plasma in such reactions^. 

One of the aims of the present paper is to show in fact that the classical and quantal 
bremsstrahlung and mini-jet sources of mid-rapidity gluons are actually equivalent up to 
form factor effects over a continuous range of ~ ^yx regime. In addition we explore 
the limits of the validity of each approximation and compute numerically the energy and 
atomic number dependence of the McLerran-Venugopalan density parameter, x- This pa- 
rameter is the total color charge squared per unit area of partons with rapidities exceeding 
some reference value. 

The calculation of this paper checks that there is a region of overlap between the 
classical and quantum computation. The quantum calculation should be valid at large 
transverse momenta. The classical calculation is valid at momenta Aqcd ^ k± ^ ^/s. 
Most of the gluons are produced in the region appropriate for the classical calculation. It 
is well known that perturbative calculations of gluon production are power law sensitive 
to an infrared cutoff. The classical computation has this infrared cutoff built into the 
calculation and may ultimately lead to a proper computation of gluon production. The 
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region where we can compare the calculations is at k± much greater than this cutoff. 

The plan of this paper is as follows: In section 2, we review the classical derivation of 
induced gluon radiation in the McLerran-Venugopalan model. We correct the treatment 
in [0 of the contact term in the classical equations of motion for a single nucleus. We 
extend further that derivation to treat properly the renormalization group corrections 
to the density parameter x- Those corrections |]TD[ increase significantly the color charge 
squared per unit area relative to the contribution of the valence quarks thusfar considered 
10, |TT], |T^. We also correct omitted factors of 2 and 27r in the original computation^]. 

In the third section, we review pQCD based derivations of induced gluon radiation[0]- 
Q. We show that the classical result agrees with the quantal results of Bertsch and 
Gunion.]^ and also with the GLR formulation]^] if DGLAP evolution of the structure 
functions is assumed. We then compare the bremsstrahlung distribution to the mini 
jet distribution and show that while the latter dominates at high transverse momentum 
p± ^ ^/x, the former dominates at p± <^ y/x- However, there is a continuous range of 
momenta ~ -^/x where both results agree at the level of 20 — 40%. 

In the fourth section, we compute the parameter x of McLerran-Venugopalan 
model. We find that due to the rapid rise of the small x gluon structure functions, y/x 
approaches on the order of 1 GeV by LHC energies for A ~ 200. Possible implications 
and further extensions of this model conclude the paper. 



2 Classical Yang-Mills Radiation 

The basic assumption of the classical approach that follows is that the coupling strength 
is small at the scale 

. 9 1 dN .9 

The parameter is the number density of gluons per unit rapidity per unit area. . 

The gluon distribution was shown in Ref. to solve an evolution equation which 
in various limits is the BFKL equation|T^, the DGLAP equation |]T6[, or its non-linear 



generalization]^. In the ultra-relativistic domain, also the rate of change of the multiplic- 
ity per unit rapidity 

is small. The smallness of this parameter means that if we compute the gluon distribution 
in a small region Ay ~ 1 around ?/ = 0, then the source of those gluons is dominated by 
hard partons with rapidities much larger than 1. These hard partons can be integrated 
out of the effective action which describes the color field source at ?/ ~ 1, and they lead 
to an effective external classical static source for the gluon field. 

Since this can be done at any reference rapidity, the classical gluon field may be 
thought of as arising from a rapidity dependent classical source. For a single nucleus 
moving near the positive light cone, we have 

D^F^'' = g^5^+p{x-,x^) (3) 
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with the source approximately independent of = (t + z)/y/2. Two types of rapidity 
variables must be differentiated. In the classical equation of motion, the coordinate space 
rapidity is relevant as defined by 

y = \nl/x~ = Vproj - \n{x'/Xp,^j) , (4) 

and x~ = {t — z)/y/2. The momentum space rapidity is, on the other hand, 

y = ^ln{p^/p-) (5) 

where 

p^ = l=(p^±p^) (6) 

are the conjugate momenta to x^. For a hadron with p+ = Pproj, we define x~j,^j = 

eXp{-yproj) ~ l/Pproj- 

The coordinate space rapidity is of the same order as the momentum space rapidity, 
since by the uncertainty principle Ax~ ~ l/p~^- Qualitatively, these rapidities may thus 
be thought of as interchangeable. On the other hand, the classical equations of motion are 
described by coordinate space variables, and we must use the coordinate space rapidity. 

In the McLerran-Venugopalan model, the source rapidity density 

p{y,x±) = x~p{x~,x±) , (7) 
is assumed to be a stochastic variable which is integrated over with a Gaussian weight, 

J[dp]exp J dyd'^x±--^TTp'^{y,x_i)^ . (8) 

This Gaussian assumption ignores correlations which we will see later are needed to regu- 
late the infrared singularities. Here fi'^{y) is the average charge squared per unit rapidity 
per unit area scaled by 1/{N^ — 1) 

Note that this p'^{y) specifies the rms fluctuations of the charge transverse density at 
a fixed rapidity. The quantity analogous to the rapidity independent p used in is 
the integrated transverse density of color charge arising from hard partons exceeding a 
reference rapidity. To emphasize this distinction we denote this quantity by 

X{y) = / dy'p'iy') . (10) 

This quantity will related below to the integrated gluon structure function. 

The solution to the above equations may be found in the light cone gauge by assuming 
that 

A± = , A' = A\y,x^) 
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The index i = 1,2 ranges over only the two-dimensional transverse coordinates. The field 
solves 

-D,-^A' = g'p{y,x^) . (11) 



Equation (|TTD is solved by lettingpf-P 



A{y,Xi_) = -[Pe''p'-°^ V Pe^p'-°j (12) 
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In this equation, P denotes path ordering along the integration in rapidity. 

If we now change variables (with unit determinant in the integration over sources) 

Fe ■'yp^oj ] p [ Pe •'yp^oj \ ^ p ^ (^13j 



then A is seen to obey the two dimensional Poisson's equation 

- ViA(2/,x^) =/p(y,x^) (14) 

Note that due to the expected slow variation of the source density as a function of 
rapidity, the field is almost constant in y. At zero rapidity, therefore, the field may be 
taken approximately as 

Ai{x-,Xi_) = 9{x~)at{x±) (15) 

where 

cti{x±) = - \Pe I Vj I Pe ^p"-"^ I (16) 

This is the non-abelian Weizsacker- Williams field of the projectile nucleus which must 
still be averaged over the ensemble (§). 

In order to generalize the above solution to the case of two colliding nuclei, we use the 
same variables as above for the projectile nucleus propagating in the +z direction. For 
the target nucleus propagating in the —z direction, we use the rapidity variable 

1/ = -|/cm + ln(X([/x+) (17) 

Here we denote the projectile rapidity with the center of mass rapidity as, ycm = Vproj- 
We will also henceforth use the index + to refer to y > and — to y < 0, when no 
confusion will arise with respect to light cone variable indices. 

In the neighborhood of y = 0, we can ignore the small rapidity dependence of the 
fields. The solution to the equations of motion in the 

x+A- + X-A+ = (18) 

gauge is approximately given by 

A^ = ±x^e{x+)e{x-)(3{T, x^) (19) 
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and 



Here r = \/t'^ — is a boost covariant time variable. (Note that the above notation 
corresponds to /3 = a and a\ = ai± of 
The fields 

af = 1 (^Pe^/L^'^^'^^^'^^^y V. (^Pe^/±.™'^^'^(^'^^)) (21) 

where 
and where 



-VlA{y,x^)=g'p{y,x^) (22) 



p{y,x±) = 9{y)p^{y,x±) +9{-y)p {y,x±) (23) 

Notice that in this solution, the fields af are two dimensional gauge transforms of 
vacuum fields. Their sum is of course not a gauge transform of vacuum fields, and therefore 
the solution cannot continue into the region x^ > 0. There is in fact a singularity in the 
solution at x"*" = and x~ = 0, at x"*" = for x~ > 0, and at x~ = for x~^ > 0. For 
x^ > 0, the form of the fields chosen above solves the classical equations of motion. In 
this region, the solution is a boundary values problem with the boundary values specified 
on the edge of the forward light cone. 

To determine these boundary values, we solve 

= g^p (24) 

and 

D^F^"' = (25) 

First we find the singularities of Eqn. |2^. In this equation, there is a S{x~^)S{x'^) 
singularity, that is a singularity at the tip of the light cone. The absence of such a 
singularity requires that 

|r=o= 0!f{x±) + a~{x±) (26) 

There are also singularities of the form 6{x^) for x^ > 0. The absence of these singularities 
requires be analytic as r ^ 0. 

The solution for the Eqn. is a little trickier since there are some potentially singular 
contact terms. It can be shown that if the fields af are properly smeared in rapidity so 
that they really solve the equations of motion in the backwards light cone, then all such 
contact terms disappear. We find that f3 must be analytic at r = and that 

P\r=o= '^[af,ar] . (27) 

The boundary conditions are precisely those of Ref. 0. They have been rederived here 
to properly account for any singularities in arising from contact terms in the equations 
of motion. These contact terms when properly regulated do not affect the boundary 
conditions. 
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We now construct an approximate solution of the equations of motion in the forward 
hght cone. We do this by expanding around the solution which is a pure two dimensional 
gauge transform of vacuum which is closest to a+ + a^. To do this, we introduce the 
projectile and target source charge per unit area at a reference rapidity y as 

r±ycm 

q^{y,x^) = ± dy'p{y',x^) (28) 



and 

, r±ycm 

ri^{y,x^)=± dy'A{y\x^) (29) 

■'y 

Note that 

{Qaiy, x±)(ibiy^ ^'±)) = x^{y)Sa,bS^{x± - x^) (30) 

in terms of x^(l/) defined as in (p^. 
By direct computation, as in 



and 



at = V,v^-'-[r]^,V,v^] (31) 

= /^g^ (32) 

The sum of a"*" + a~ can be written as a pure two dimensional gauge transform of 
vacuum plus a correction as 

af + a~ = + 5q;° (33) 

where 



a. 
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and where 

Sa^='-{[v-,^^V-'] + [v^,^^V-]} (35) 

This decomposition into a gauge transform of the vacuum is accurate up to and including 
order g"^. 

Now we expand af = + 5af(r, x±). Both 6a^ and /5 are the small fluctuation fields 
corresponding to radiation. We find that and /3 solve exactly the same equations as 
were incorrectly derived in Ref. P|. So even though the original derivation was incorrect, 
the final result remains fortunately valid. 

In Eqn. 42 0, a factor of 2tt was however omitted, and as well in Eqns. 45, 47, 49 
and 50. In addition, in going from the first of Eqns. 49 to the second, a factor of 1/2 
from the trace was omitted. 

The final result corrected for the above factors and generalized to include the source 
of hard gluons is 

dN ^2^g'x^iy)x-iy) n,{n! - i) f d\^ 

— ttK 



dyd^k^ (27r)3 kl J (27r)2 g2 (q^ _ k^)2 

:^y) 

(2vr) 



ttR^ ^ / f ,4 -L{k^,X) . (36) 
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The q± = and q_L = divergences arise in the above classical derivation because of 
the neglect of correlations in the sources ensemble. A finite logarithmic factor , L{k±, A), 
is obtained only if we include a finite color neutralization correlation scale, A. 

This scale arises from dynamical screening effects and may be seen in models such as 
the onium valence quark model of Kovchegov|ll] as developed in [12]. In the classical 
calculation, this cutoff appears after averaging over various values of the background 
charge density. |jlO| The cutoff scale turns out to be A ~ ay/x- Below this cutoff scale, the 
factors of l/fc^ moderate and become of order ln{k±). This cutoff scale acts somewhat like 
a Debye mass, although this is not quite the case since the logarithmic dependence implies 
power law fall off in coordinate space whereas a Debye mass corresponds to exponential 
decay. In any case, for evaluating L{k±, X) at k± » X the precise form of the cutoff 
is unimportant, only that the l/k'j_ singularities in the integrand are tempered at some 
scale. This is because logarithmically divergent integrals are insensitive in leading order 
to the precise form of the cutoff. The generic form of the logarithmic factor is therefore 
expected to be of the form 



L{k±,X,y) 



2n 



q± - kj 



(37) 



where JF is a suitable form factor. In a dipole form factor was considered. A gauge 
invariant screening mass was considered in . Such dipole form factors lead to 



L{k^,X) = klj 



27r 



log(A;l/A^) 



(38) 



(gi + A2)((qx-k^)2 + A2) 

where the logarithmic form is remarkably accurate for k±/X > 2. A finite but nonlogarith- 
mic form of L can also arise if other functional forms for the form factors are considered 
as in 
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It is also important to stress that in any case, the above classical derivation neglected 
nonlinearities that can be expected to distort strongly the above perturbative solution 
in the k\ < o?x region. Thus, the solution should not be extended below A ~ in 
any case. In future studies, it will be important to investigate just how the full nonlinear 
Yang-Mills equations regulate these infrared divergences. 



2.1 Classical Color Current Fluctuations 

For two colliding nuclei the effective classical source current for mid-rapidity gluons is 
assumed to be 

j,^(x) =5'^+5(a;-)g+(0,x^) + 5'^-5(x+)g,-(0,xx) (39) 

where (g^) = but the ensemble averaged squared color charge density of each of the 
A^^ — 1 components is given by x^(0) as in (^). 

In ref . [pT] , , x was estimated using the valence quark density and with the classical 
color density interpreted as a color transition density associated with the radiation of a 
color a gluon 

N 

g'^(x^) = ^(T„"),,,52(xx-x^„) (40) 

n=\ 
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where the sum is over the valence quarks, and T° is a generator of dimension dn ap- 
propriate for parton n. In this interpretation, products of color densities involve matrix 
multiplication and the ensemble average leads to a trace associated with averaging over 
all initial colors of the valence partons and a summing over all final colors. Therefore 

k"(x±)) = (41) 
since TrT'^ = in any representation while 

(g<^(x^)g^(x'J) = -rTr(Ty:)n{-^^)5\^^ - x/) (42) 

n=l ""n. 

where n(x_L) = (5(x_l — x^(n))) is the transverse density of partons of type n. From now 
on we assume identical projectile and target combinations and fix ?/ = so that we can 
drop the distinction between ± sources and the rapidity variable. 

Taking into account both the valence quark and hard gluon contributions in the nuclear 
cylinder approximation used in p|, the relevant x = X^(0) parameter is therefore given 
by 



c 



where the transverse density of quarks is ng(x^) = Nq/TrR"^ and the gluon transverse 
density is ^^(x^) = Ng/nR'^, 

Because this interpretation allows for complex color (transition) densities that do 
not arise in the classical limit, it is useful to show that it can also be derived from a 
more conventional classical Yang-Mills treatment. For that purpose we use the Wong 
formulation of classical YM kinetic theory ||T5|. In that formulation, the parton phase 
space is enlarged to incorporate a classical charge vector A"(r) in addition to the usual 
{x'^ {t) , (t) = mu^ij)) phase space coordinates. The phase space density, /, obeys the 
Liouville equation 

'^-/(x(r),p(r),A(r)) = (44) 



dr" 



with dx^/dr = and 



dii^ 

m— = gu^FrA' 

^ = gP^^u.A'^^A^ = -iir^e^) A". (45) 
dr \ / ac 

where {T^)ac = if"-^'^ are the generators in the adjoint representation and 9^ = gu^A^^. 
The color current gi^{x) in this kinetic theory is computed via 

f'^x) = J dTu>'{T)A^{T)6\x - x{t)) (46) 

The color charge vector precesses around the local field but its magnitude remains 
constant. Its length is fixed by the specified color Casimir C2 = J2aAl- In the ultra- 
relativistic case with p^/po ~ 1, the current reduces to eq . (p9|) with the transverse density 

g'^(x^) = Y: U,cA:{t,)6\^^ - x^J (47) 
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where the unitary U = P exp{—ig Jq dsT^Ufj_A^^^{x{s))} accounts for the color preces- 
sion along the parton trajectory. The ensemble average in this formulation involves an 
integration over the initial colors A5^(ro) with a measure 

dAn OC n dA'J{AX - C2n) (48) 
c=l 

normalized such that / dA^ = 1 and thus 

I dAr^AlAi = S'^'C^n/dA (49) 

Because U is unitary, this leads to the same expression for the color charge squared 
correlation parameter x as eq.(^3|). 



2.2 Yang-Mills Radiation Distribution 

Inserting the above expression for x iiito the classical formula for radiation, we obtain 
dN _ I ( N, ^ N.NgY 2/ N,{N^-1) f d\± 



dyd'^k± TiR^ \2N^ N^-lJ (2vr)3 kl J (27r)2 g2 (q^ _ k 



\2 



^ (CpN, + CAN,)f^'^^^L{k^,X) (50) 



7ri?2 9JJ (27r)4 kl 

If only valence quarks are included then this reduces to 

In the opposite limit, if only hard glue is included, the radiation distribution reduces to 

dN f2g^NA ( C\\ 1 



dyd^k^ 7ri?2 I {2txY / \dA , k\ 



L{k^,\) (52) 



Note that the color factor in the second brackets marked el is that associated with the 
elastic scattering of two partons 

C^m = ( j-JrT:T!:) (^TrT^T^) = , (53) 



so that C^' = 2/9,9/8 for qq,gg. The elastic Rutherford cross section is in this approxi- 
mation 

~ (27r)2 ] ~ Tt'' dA ■ ^ ^ 

The infrared divergence is regulated by the color screening scale A or form factors as in 
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The geometrical Glauber factor in Eqns. ( counts the average number of binary 

parton-parton collisions per unit area in 6 = collisions of cylindrical nuclei. More 
generally, 

Tab{^) = / d^XPn/A{^) / dZB 0-n.mPm/B{^± " b, Zfi) (55) 

For 6 = collisions of cylindrical nuclei this reduces to 

N N 

T""^(0) = (56) 



Therefore we can write 



T'""iO)'!^ (57) 



where 



dyd'^k^ dyd^k^ 



dyd^k^ ""^ V (27r)3 ) klJ (27r)2 (q^ - kj 

3 Quantum Radiation 
3.1 pQCD Bremsstrahlung 



We compare (|58| ) with the quantum radiation formula derived in 0. In the = 
gauge and for gluon kinematics k = [xP^, k'j_/xP^, k^] with x <^ 1, the three dominant 
diagrams sum in the small momentum transfer limit to 

zM{nm ^ g{k,e,c)) = [TIT:]T:^ {^-^^ [2ge^ • |^ - (59) 

Taking the square and averaging over initial and summing over final colors, one finds that 

da I C^^ATTa^\ ( aNc q\ \ 



dq^dyd^k^^ \ y y tt^ k'j_(k± — q_|_' 



2 



dq'^ dyd'^k_L 

This is the basic factorized form of the soft QCD radiation associated with elastic scat- 
tering. Integrating over the elastic momentum transfer q± yields 

da f d'^q± ( C^^Ana^ \ ( aNc q]_ \ 



dyd'^k± J \ q± / \ k'j_(k± — q±) 

, 2g'Nc 1 f d\^ 1 
""^ (27r)3 kl J (27r)2 ql{kj_ - q^^2 ) 

This is exactly the same as the classical result in (^). 



11 



In Ref. 0] the nir —* g cross section was computed taking a dipole form factor into 
account with the result 

da-^-s (Caoc'\ /•rf^q±22F^(gi)F^((qx-kx)2) 

(62) 



dydkl \7i^klJ J (27r)2 ql{q^ 
where 

4g 



P-(l) = . 2I 2 (63) 

Again we can read off the elementary qq ^ g cross section by dividing by the number of 
parton pairs = 4 in this reaction and neglecting interference by setting = 1. This 
leads to 

daii^a 1 f2g^NA 1 f d^q^ 1 



dydkl 4 V (27r)3 J kl J (27r)2 ql{q^ - k 



(64) 



where the first factor 1/4 is just the large N^. limit of C^^ — > 1/4 used implicitly in eq.(17) 

of i- 

3.2 Comparison with GLR formula 

It is also of interest to connect the classical YM formula with the pp ^ g formula of 
Gribov, Levin, and Ryskin (GLR)mand used recently in ref.Qto compute mid-rapidity 
gluon production at LHC energies: 

da aNc /• 2 fixi,ql)f{x2,ik±-q±f) 

- Knotty / " q± jTi v? ' (65) 



dyd'^kj_ 7r2/c]_ J 9i(k± — qi_y 

where 



and 



f{x,Q') = —^^xG{x,Q') (66) 
d log Q-' 



xi ^ X2 ^ xj_ = k^l\fs (67) 



are fractional momenta which are assumed to be small. In this relation, the radiation 
resulting from the fusion of two off-shell ~ 2/2 ~ gluons is estimated. Unfortunately, 
there is variation in the literature as to the magnitude of the factor |Q. This is partly 
due to variations in the definition of /(x, Q'^\ We find below that in order to reproduce 
the perturbative QCD and classical Yang-Mills result, eq.(^,^), we must take 



= ^ 5 . (68) 

From private communication with E. Levin, this factor is required if / is defined as in 
ref.P] via (|66|). This implies that the results quoted for the BFKL contribution to minijets 
in taking Kn = 1 are approximately a factor 5 too small. With the value in eq. (|68D , 
the BFKL and conventional mini-jet rates would coincide more closely. In section 3.3 we 
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argue that at least in the asymptotic domain these distributions should in fact coincide 
over a range of ~ a/x- 

To compare to the GLR formula with the classical result, we approximate the 
evolution using DGLAP evolution p6[ 

In the small x semi-classical domain 

P,g{x) ^ 2Njx (70) 
Therefore, we have the approximate relation at high 

fix, Q2) ^ ^ dx' G{x', Q2) = ^Ng{x, Q2) (71) 

Consequently, from (|66| , |7TD 

da aNc r ^ d xG I d xG \ 



dyd'^'k± 7r'^kj_ J dq\ \ dq 



-,12 



(kx-qx)2 



aN, a^N^ r iV,(xi, gi)iV,(x2, (k^ - q±) 
^N _^,o _o / d q±- 



(72) 



Equation (^) reduces to the classical expression ( |5^ ) if we approximate the integral by 
factoring out the integrated gluon numbers at the average scale, ~ k^, divide by nR'^, 
and take the normalization factor Kj^ from (BSf). 



3.3 Matching 2 ^ 3 to 2 ^ 2 

Up to this point, we have shown that the classical and quantum bremsstrahlung formulas 
agree for the 2^3 process up to specific form factors. The problem addressed in this 
section is the relationship between the bremsstrahlung spectrum and the mini-jet spec- 
trum based on the pQCD factorized 2 — > 2 processes. Recall|P] the factorized differential 
cross section for two gluon jet production with transverse momenta , ±k_L, and rapidities 
yi and y2, is given by 

- — - — — — = K xiGA{xi,k\)x2GB{x2,k']_) — , (73) 

"2/l"y2" k_L TT dt 

where Xi = x±{exp{yi) + exp{y2)) and X2 = x±{exp{-yi) + exp{-y2)) , with x± = k^/^s, 
and where the pQCD gg — » gg cross section for scattering with t = —kj_{l + exp(t/2 ~yi)) 
and 2/2 — 2/1 = 2/ is given by 

^ ,W(l + e^+e25^ 
dt ki {evI^ + e-yl^Y ' ^ ' 
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This reduces to the naive Rutherford expression eq. only if the unobserved gluon has 
a rapidity \y\ ^ 3. For \y\ < 1, the exact form (ff^) is 27/64 0.42 smaller than the 
Rutherford approximation. 

We concentrate here only on the dominant gluon-gluon contribution for symmetric 
systems, A + A, with G = Ga- The inclusive gluon jet production cross section is obtained 
by integrating over y2 with yi = y and k± fixed. For an observed midrapidity gluon with 
?/ = 0, —y* < y2 < y*, where exp(— y*) = x±/{l — x±), we must evaluate 

I{x±,k]_) = dy2 XiG{xi, k]_)x2G{x2, k\) ^^^^^^ ^ £-^2/2)6 ' 

(75) 

In the Rutherford approximation, implicit in the classical approximation, we neglect the 
y dependence of (^) and therefore approximate I by 

/J/* 
dy2 XiG{xi,kl)x2G{x2,kl) (76) 
-y* 

^ 2x±G{x±,k'i) / dx2G{x2,k']_) 
^ 2x^G{x±,kl)Ng{x±,kl) . 



For xG oc x^ (l—x)'^ with 6 ~ 0.2, 7 ~ 8.5, as HERA dataP7|, ^ indicate in the moderate 
~ 5 GeV^ range, the last approximation to Ir is found to agree remarkably within 
10% of the numerical integral of the first line as long as Xj_ ^ 0.01. However, for k'j_ = 5 
GeV^, the neglect of the rapidity dependence of dagg/dt in the Rutherford approximation 
leads Ir to overestimate / by ~ 55% at RHIC energies {x± ~ 0.01) and by ~ 34% even at 
LHC energies x± ~ 0.001. This is due to the factor ~ 2 suppression of the pQCD rate in 
the — 1/2I < 1 range. On the other hand, next to leading order corrections modify ([73|) 
by a factor iiT ~ 2 in any case, and the next to leading order corrections to the classical 
formula are not yet known. Since neither the mini-jet nor the classical radiation can be 
determined at present to better than ~ 50% accuracy, the following simplified Rutherford 
formula for the single inclusive pQCD mini-jet distribution is adequate: 

2Ngix,t)xGix,t)(^] (77) 



R 



dydt \ dt 

or 

dN_ ^ 2iV,(xx,t) ^^^^^^^ fda^\ ^^g^ 

dydt ttR"^ ' \ dt ) 

In order to compare the above mini-jet distribution with the classical bremsstrahlung 
result (0), we need to replace the factor in ([5^) by Ng{x±, q'j_)Ng{x±, {k± — q±Y) and 
move that factor inside the logarithmic integrand. This generalization is essential since 
the effective classical source due to hard gluons depends on the x± and scale of resolution 
of the radiated gluon. This requires that k'j_ be sufficiently large so that the variation of 
the structure function with that scale be small. In this case, the classical bremsstrahlung 
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formula generalizes into the GLR form (^) 



dydt ttR^ " dt (27r)2 ql (k^ - qx)^ 

1 ,o.dat r^l , o d 



dC Z*^ I n 

- 2iV(x^, fcl)-f / " gl) , (79) 

Trit^ at JO agj^ 

where in the last step we used the DGLAP evolution (0). Thus, we recover the same 
minijet formula as (|78|). 

The use of the DGLAP evolution is essential to prove the duality between classical 
bremsstrahlung and the conventional minijet distributions. We note that in order for 
corrections to the classical result to remain small, it is necessary that a{x) ^(fc^/a^x) 
1. Recall that in the classical analysis, a is always to be evaluated at some scale of order 
^/x ^ ^QCD- This requirement is therefore that k± ~ y/x- If this is satisfied, then the 
formulae should agree in the <^ 1 regime. 

We see therefore that all the formulae used for hard scattering agree with the classical 
result in the range of momentum a^x -C /c^. ^ X- This range of momenta is outside 
the typical scale kl ~ a^x which non-trivial behaviour of the transverse momentum 
distributions is expected on account of screening. In the region of smaller k±, the full 
nonlinearity of the Yang Mills equations must be taken into account. At large k_i_ > y/x, 
the hard scattering pQCD formula properly sums up higher order DGLAP corrections to 
the classical formula. It is important that there is a range of momenta where the classical 
and hard scattering results match at the level of ~ 50%. 

4 Estimate of s, Q^) 

We turn finally to the estimate of the McLerran-Venugopalan scale density x in the range 
of A and s in future RHIC and LHC experiments. 

4.1 Valence Quark Contribution 

The initial assumption in [|l| and further developed in [11, 12] was that for A ^ 1, 



the valence quarks could provide a very high density of hard color source partons for 
which recoil effects are negligible and thus treated classically. In the nuclear cylinder 
approximation, the transverse density of valence quarks is simply 

N^A N,A^/^ , , 

where ro = 1.18 fm. Since each quark contributes with a color factor Cp/dA = l/2Nc the 
valence quark contribution to the color charge squared density 

Xvai = 7r-2 = 0-07 GeYf < AlcD , (81) 
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where the bound arises because only A < 200 beams will be available. One would need 
astronomical A ~ 10^ to reach ^^Xvai = 1 GeV because of the extremely slow A^^^ 
growth. Thus, the valence quark contribution is in practice too dilute to contribute into 
the perturbative domain. Fortunately, the non-Abelian Weizsacker- Williams field has a 
very large number of "semi-hard" gluons with x ^ x± in the s — >■ oo limit. The question 
then is how large does s have to get in order to push x into the perturbative regime. 



4.2 The Hard Gluon Source 

The number of hard gluons that can act as a classical source for midrapidity gluons 
depends not only on x± but also the resolution scale Q 

Ngix±,Q^) = C G{x,Q^)dx . (82) 

Each gluon contributes to x with a color factor CaI^a = 3/8. In order to get an 
upper bound, we will neglect possible nuclear glue shadowing and assume that Ng oc A. 
There could be some suppression of the low x gluon number in nuclei due to shadowing as 
observed for nuclear quarks. In the McLerran-Venugopalan model, only a mild logarithmic 
~ ln(y4) modification of A^g oc A is expected. Including valence and sea quarks and 
antiquarks as well as gluons leads then to 



TT Tq Jxo V6 8 / 



The lower bound, Xq, is determined up to a factor of ~ 2 by the minimum momentum 
fraction needed to justify the neglect of recoil associated with the radiation of a midra- 
pidity gluon with a given k^. In the estimate below we vary that bound between xq = x^ 
and 2x^. For xg{x) oc Xjx^ with b ~ 0.2 — 0.3, this leads to an uncertainty, bxlx ~ 
well within the present overall normalization uncertainties of the small x gluon structure 
function. 



In order to compute x we use the GRV95 NL0(M5') parametrization|p!9| of the nu- 
cleon structure functions. Figure 1 shows how the parametrization of the gluon structure 
compares to preliminary "data" at = 7 GeV^ from HERA P^, obtained via a 
DGLAP analysis of the scaling violations from F2(x,(5^). Also shown is the BFKL-like 
parametrization of the gluon structure used in for comparison. Both the GRV95 and 
the BFKL parametrizations significantly overestimate the moderate data of interest 
here at x < 10~^. The preliminary = 20 data from II1]TB| (not shown Fig.l) also lie 



below the GRV95 parametrization. For our purposes, it is only important that the use of 
GRV95 and the neglect of gluon shadowing should lead to a reasonable upper bound on 
X- 

As discussed in the previous sections, the classical regime extends up to k± ^ Qym 
where 

1/6 /]^Q-4\^/2 



QyaM, s) = VxiA, s, Ql,,) ~ (1 GeV) (^— j [-^ 



(84) 
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and 6 ~ 0.2 — 0.3. In principle, this must be determined self-consistently given the scale 
dependence of the glue. In practice, as shown below, Qym is only weakly dependent of 
the reference scale if its above ~ 2 GeV^. The approximate formula for Qym summarizes 
the numerical results below. 

In Figure 2, Qym is plotted for y/s = 0.2, 6.5, 100 ATeV for heavy nuclear beams 
with A = 200 as a function of the scale Q with which the GRV95 structure functions are 
evaluated in (0). The upper solid curves for each energy correspond to (^) with xq = Xj_. 
The lower curves are obtained by increasing the lower cutoff from x± to 2x±. The two 
long dashed curves at the bottom show the contribution of only the valence quarks to 
Qym at y/s = 0.2, 100 ATeV using xq = x±. The curves show that the hard nuclear 
glue dominates for finite nuclei at all collider energies. Note also that Qym is remarkably 
independent of the reference Q scale because of a compensation of two competing effects. 
The increase of xG with Q is compensated by its decrease with increasing value of the 
minimum hard fraction xq ^ Q/^ys contributing to the classical source. 

At RHIC energies, the boundary of the classical regime remains rather low ( ^ 500 
MeV) because the relevant x range, x± > 0.005, is not very small. By LHC energies, 
on the other hand, gluons down to x± ~ 0.0001 can contribute, and the classical Yang- 
Mills scale increases to Qym ~ 1 GeV. Note that to double the Qym scale at a fixed 
energy would require an increase of A to 2^ x 200 if shadowing can be ignored. To double 
the value of Qym at fixed A requires decreasing x± by a factor 2~^° ~ 10~^. Although 
asymptotically the scale of Qym becomes arbitrarily large, this asymptotic behaviour is 
approached slowly. We conclude that at RHIC energies the classical Yang-Mills radiation 
dynamics is likely to modify mainly the nonlinear, nonperturbative beam jet regime. In 
that regime the perturbative analysis must certainly be extended into the full nonlinear 
regime via detailed numerical simulations. By LHC energies it appears that the classical 
Yang-Mills radiation begins to overlap into the perturbative minijet domain with fc^ ~ 1 
GeV. 

The very small x and very large A limits, where perturbative classical radiation can 
be computed, provide an novel calculable theoretical limit. It provides qualitatively use- 
ful insight at RHIC energies and may be semi-quantitative already at LHC energies. In 
future studies it will be especially important to extend work with this model into the non- 
linear regime to clarify the mechanisms for color screening in A -|- A reactions at the lower 
k± ~ aQYM scale. Present estimates for initial conditions in A -|- A based on mini-jet 
pQCD analysis^, ^ |] vary considerably because of the necessity to introduce a cutoff 
scale, po ~ 1 ~ 2 GeV, to regulate the naive infrared divergent Rutherford rates. That 
cutoff has thus far been estimated either (1) phenomenologically by imposing observed 
constraints from extensive pp, pp vr, K, p, X systematics as in 0, |] or (2) using ki- 
netic theory estimates which are sensitive to formation physics effects. One of the 
great theoretical advantages of the classical Yang-Mills approach is that the long wave- 
length nonlinear dynamics involving pre-asymptotic field configurations can be taken into 
account (at least numerically) without invoking kinetic theory or formation physics as- 
sumptions. In the theoretical uQym ^ 1 GeV domain, that physics may be accessible 
using perturbative techniques. In the experimentally accessible uQym < 1 GeV regime, 
numerical solutions of the Yang- Mills equations, as for example in [^, are likely to pro- 
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vide additional insight into that problem. The classical Yang-Mills model ^ is one of 
the practical tools at present to approach the study of asymptotically high energy reac- 
tions, where many unsolved and interesting theoretical problems remain. 
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5 Figure Captions 

Figure 1: The GRV95 NLO[||] and BFKL-likeg] parametrizations of the gluon 
structure function, xG{x,Q'^), for Q'^ = 7,20 GeV^ are compared to preliminary ZEUS 



"data" 13 from HERA. 



Figure 2: The classical Yang-Mills scale, Qym, from (|83|) is shown for A = 200 nuclei at 
collider energies y^s = 0.2, 6.5, 100 ATeV as a function of the reference scale, Q, used to 
evaluate the GRV95|p!9| structure functions. Upper curves and lower curves for each 



energy correspond to taking the lower cutoff scale xo = x±, 2x±, respectively. The 
bottom two dashed curves give the valence quark contributions at = 0.2, 100 ATeV. 
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GRV95 vs BFKL Parametrizations of xG(x,Q ) 
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Figure 1: The GRV95 NLO[T^ and BFKL-like|^ parametrizations of the gluon structure 
function, xG{x,Q'^), for Q'^ = 7,20 GeV^ are compared to preliminary ZEUS "data" [p!^ 
from HERA. 



20 



Classical Yang-Mills Scale for A=200 

Qym = A'^' (0.1 GeV) Mx (q(x,Q')/6+3g(x,QV8) 
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Figure 2: The classical Yang-Mills scale, Qym, from (^) is shown for A = 200 nuclei 
at collider energies y/s = 0.2,6.5, 100 ATeV as a function of the reference scale, Q, used 
to evaluate the GRV95|T^ structure functions. Upper curves and lower curves for each 
energy correspond to taking the lower cutoff scale Xq = x±, 2x±, respectively. The bottom 
two dashed curves give the valence quark contributions at y/s = 0.2, 100 ATeV.. 
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